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ABSTRACT 

We investigate the effects of finite sky coverage on the spectral resolution A£ in the 
estimation of the CMBR angular power spectrum C e . A method is developed for ob- 
taining quasi-independent estimates of the power spectrum, and the cosmic/sample 
variance of these estimates is calculated. The effect of instrumental noise is also consid- 
ered for prototype interferometer and single-dish experiments. By proposing a statistic 
for the detection of secondary (Doppler) peaks in the CMBR power spectrum, we then 
compute the significance level at which such peaks may be detected for a large range 
of model CMBR experiments. In particular, we investigate experimental design fea- 
tures required to distinguish between competing cosmological theories, such as cosmic 
strings and inflation, by establishing whether or not secondary peaks are present in 
the CMBR power spectrum. 

Key words: cosmology: miscellaneous - cosmic microwave background 



1 INTRODUCTION 



The cosmic microwave background radiation (CMBR) is 
one of the most promising links between astronomical ob- 
servation and cosmological theory. A significant amount of 
experimental data already exists and further observational 
progress is expected in the near future. The experimental 
success in this field has prompted a large amount of theo- 
retical effort on two closely related fronts. Firstly, theorists 
strive to assess the impact current observations have on the- 
ories of the early Universe. Secondly, the design of future ex- 
periments is guided by what are believed to be crucial tests 
of our theoretical prejudices. 

A particularly active field of research in CMBR physics 
are the so-called Doppler peaks (Hu & Sugiyama 1995a,b). 
These consist of a series of oscillations in the angular power 
spectrum of CMBR fluctuations C l predicted for most infla- 
tionary models. They are predicted in the multipole range 
100 & i Ss 1500, corresponding to angular scales 0.05 ^ 
6 ^ 1 degrees. Experimental measurement of the Doppler 
peaks' positions and heights would fix at least some combi- 
nations of cosmological parameters (e.g. Ho, f2o etc.) which 
are left free in inflationary models ( Jungmann 1995) . Future 
experimental projects are usually designed with these goals 
in mind. In particular, models with flo = 1 are of special 
interest. 

In this paper we shall consider another theoretical con- 
text. For a long time inflation (Steinhardt 1995), and topo- 
logical defects (Kibble 1976; Vilenkin & Shellard 1994), 



have stood as conceptually opposing alternative scenarios for 
structure formation in the early Universe. It was shown by 
Albrecht et al (1996) and Magueijo et al (1996), that qualita- 
tive aspects of Doppler peaks should reflect this conceptual 
opposition. In particular the absence of secondary Doppler 
peaks was proved to be a robust prediction for a large section 
of defect theories. Interestingly enough this statement relies 
only on the role played by causality and randomness in these 
theories. Cosmic strings were unambiguosly shown to fall in 
this category. However, this may or may not be the case for 
textures (Crittenden & Turok 1995; and Durrer et al 1995). 
Therefore it appears that determining whether or not there 
are secondary Doppler peaks in the CMBR power spectrum 
offers an important alternative motivation for experimental 
design and data analysis. Answering this question is a far 
less demanding experimental challenge, which nevertheless 
will have a dramatic impact on our understanding of the 
Universe (Magueijo & Hobson 1996, Albrecht & Wandelt 
1996). 

We address this issue by proposing in Section |B| a statis- 
tic for detecting secondary oscillations, and studying how it 
fares for signals coming from various models, when mea- 
sured by using different experimental strategies. The result 
is encoded in a detection function E telling us to within 
how many sigmas we can claim a detection of secondary 
oscillations, given a particular model and experiment. We 
consider signals coming from standard CDM (sCDM), and 
an open CDM model (henceforth called stCDM) which is 
tuned to confuse inflation and cosmic strings in all but the 
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issue of secondary oscillations. In Sees. and ^| we set 

up a framework for computing errors in power spectra es- 
timates. We consider errors resulting from spectral resolu- 
tion limitations due to finite sky coverage, cosmic/sample 
variance, instrumental noise and foreground subtraction. We 
consider a large parameter space of experiments including 
single-dish experiments and interferometers. For single-dish 
experiments we allow the beam size, sky coverage, and de- 
tector noise to vary. For interferometers we take as free pa- 
rameters the primary beam, number of fields, and detector 
noise. We then use this framework to compute the detec- 
tion function E for sCDM and stCDM in this large class of 
experiments. 

The results obtained are given in Sections [5] and |(| and 
provide experimental guidance in two different ways. Firstly 
they allow the choice of an ideal scanning strategy (choice of 
resolution and sky coverage) given a constraint such as finite 
funding (albeit disguised in a more mathematical form) . Sec- 
ondly, one may compute the expected value of the detection, 
assuming ideal scanning, as a function of available funding. 
This provides lower bounds on experimental conditions for 
a meaningful detection as well as an estimate of how fast 
detections will improve thereafter. We summarize the main 
results in the Section 



2 OBSERVATIONS OF SMALL FIELDS 

CMBR temperature fluctuations are usually described as an 
expansion in spherical harmonics 



AT (x) \- \- 

1 = m=-f 

from which we define the CMBR power spectrum C e = 
(|a^m| 2 ). Note that in this paper we write ensemble aver- 
age quantities with the indices upstairs, and denote observed 
values of such variables in any given experiment with the in- 
dices downstairs. 

If one is looking only at small patches of the sky, the 
spherical harmonic analysis is awkward to apply, and it is 
more convenient to use Fourier analysis. For that purpose 
we perform a stereographic projection of the celestial sphere 
onto a tangent plane at the centre of our small patch, so 
that circles with colatitude 9 are mapped onto circles of 
radius r = 2tan(#/2) on the plane. We can then describe 
the CMBR fluctations AT(x) /T on this plane by its Fourier 
transform. We use the conventions 



AT(x) 
T 



and 



a(k) 



d 2 k ,,, ik . x 
2n 1 ' 



d 2 xAT(x) _ ih . t 

e 

2tt T 



If we denote the Fourier transform of the autocorrela- 
tion function of the AT(x)/T field by C(fc)/(2ir) then the 
covariance of the Fourier modes a(k) is given by 



(a(k)a*(k')) = C(k)S(k - k'), 



(1) 



where the angle brackets denote ensemble averages, * de- 
notes complex conjugation and k = |fe|. We note that, due 



to the requirement of rotational invariance, C(k) = C(k). 
From (|l|), the raw modes a(k) are therefore independent 
Gaussian random variables with variance C(k). In calcu- 
lations concerning small patches of the sky C(k) can be 
obtained by interpolating the C e coefficients in the spher- 
ical harmonic expansion with |fe| — £ (Bond & Efstathiou 
1987). An improvement can be obtained with the prescrip- 
tion C{k) = C[l{k))dk/dl with k(£) =ncot(n/£). 

We wish to study the statistical properties of the a(k) 
modes as seen by observations of a small field with a given 
observing beam. Let us describe the field by a window W(x) 
and the observing beam by B(x). The sampled temperature 
map for a single-dish experiment is then 



Ar s (a:) 
T 



AT(x) 
T 



*B(x) 



W(x), 



where * denotes convolution. The sampled Fourier modes 
a s (k) are therefore given by 



a s (k) = [a{k)B(k)]*W(k) 
d 2 k'~ 



2tt 



W(k-k')a(k')B(k'), 



(2) 



where W(k) and B(k) denote the Fourier transforms of the 
window and observing beam respectively. 

For interferometers, however, we must make a slight 
modification since in this case the Fourier domain is sampled 
directly. An interferometer samples the Fourier transform 
of the product of the sky temperature fluctuations and the 
primary beam of the antennas W(x) (which corresponds to 
a window). The positions of these samples (or visibilities) 
in the Fourier domain (or ra-plane) are determined by the 
physical positions of the antennas and the direction to the 
patch of sky being observed, and lie on a series of curves (or 
M?;-tracks). If we denote these curves by the function B(k), 
which equals unity where the Fourier domain in sampled 
and equals zero elsewhere, then the sampled Fourier modes 
for an interferometer are given by 



a s {k) = [a{k)*W(k)]B(k) 
d 2 k'~- 



B(k) 



2n 



W(k-k')a(k'). 



(3) 



The inverse Fourier transform of B(k) is the synthesised 
beam of the interferometer B(x). The a s (k) modes are of- 
ten inverse Fourier transformed to make a sky map. This 
sampled sky temperature distribution for an interferometer 



AT s (x) 



AT(x) 



W{x) 



*B(x), 



which is the convolution of the synthesised beam with the 
product of the sky and the primary beam (or window). 

From (|l|), (^) and (^|) one can then derive the covariance 
matrix of the sampled modes for each experiment (Hobson, 
Lasenby & Jones 1995) . For a single-dish observation we find 

(a s (k)a*(k')} = 



d 2 k h 
(2^) : 



W(k - k")W*(k' - k")C(k")\B{k h 



(4) 



whereas for an interferometer 
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{a s (k)a* s (k')) 



\B(k)f 



d 2 k" 
(2tt) : 



W(k - k")W*(k' - k")C(k"). 



(5) 



Ignoring, for the moment, the effects of the observing beam 
B(x) in each case, we see that finite sky coverage [described 
by the window VK(cc)] renders the observed a s (k) dependent 
random variables. Roughly speaking, if L is the typical size 
of the field then a correlation length £ ~ 1/L is introduced 
in the fe-plane. 



2.1 The sampled power spectrum 

The two-dimensional sampled power spectrum may be de- 
fined in terms of the marginal variances of the sampled 
modes as 

C s (fc) = (a s (k)a* s (k)) 

From (Q) and (^) we see that C s (fe) is given for single-dish 
experiments by 



C(k) 



(2^) 2 
= C{k)\B{k)\ 2 

and for interferometers by 
d 2 k 



d k ' \W(k-k')\ 2 C(k')\B(k')\ 2 
*\W(k)\ 2 , 



(6) 



C'(fc) 



|B(fc)|- 
|B(fe)| 2 



(M 2 

C(k)*\W(k)\ 2 



W(k-k')\ 2 C(k') 



(7) 



Ignoring, for the moment, the effects of the observing beam 
B{x) [by setting its Fourier transform B{k) to unity in (^|) 
and (pi)], we see that in both cases the two-dimensional sam- 
pled power spectrum is the convolution of the underlying 
CMBR power spectrum C(k) with |VK(fc)| 2 . We also note 
that in general C s (k) is not circularly symmetric. We there- 
fore define the one-dimensional sampled power spectrum as 
the azimuthal average C(k) = {C s (k,8))g, where is the 
azimuthal angle in the Fourier domain. 

In order to provide examples of the effects of the win- 
dow on the sampled power spectrum, we now consider two 
simple windows: a square window and a Gaussian window. 
A square window is a toy model for single-dish experiments. 
If the square window has sides of length L, then its Fourier 
transform is given by 



W{k) 



EL 

2tt 



Jo 



(¥)*(¥) 



(8) 



where jo(x) = sinx/x is the spherical Bessel function of 
order zero. A Gaussian window normalized to one at its 
peak is a toy model for an inteferometer primary beam, and 
is given by 



W(x) =exp [ -^-2- 



which has the Fourier transform 

t.2 2 

K 1( , 



W(k) 



, exp 
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Figure 1. (Top) The sampled power spectrum C s (fc,0) (dashed 
line) for a single beam experiment with a square window of size 
L = 1.8 degrees and resolution fcj, = 1000, compared to the un- 
derlying sCDM power spectrum (solid line) as observed by an 
all sky experiment with the same resolution. (Bottom) The same 
as above but for the azimuthal average sampled power spectrum 

c(k) = (C(k,e)) g . 



We note that both the square and Gaussian windows are 
real, even functions, and hence so are there Fourier trans- 
forms. 

Using these windows we may now calculate the sam- 
pled power spectrum for various experiments (assuming for 
the moment that B(k) — 1), given an underlying CMBR 
power spectrum. We choose the underlying spectrum to be 
that predicted by the standard inflation/CDM scenario with 
O = 1, ho = 0.5, and Q b = 0.05 (which we shall call sCDM). 
In order to accentuate the effects of the windows we first 
consider small fields. 

As a model single-dish experiment we consider a square 
field of size L = 1.8 degrees observed with a Gaussian beam 
with a w — 0.06 degrees (or FWHM equal to 0.14 degrees), 
so that 



B(k) =exp ( 



(9) 



where fej = 1000. For a square window the two-dimensionsal 
sampled power spectrum is not circularly symmetric, and so 
C s (fc, 8) is a function of 9. If we fix 6 — 8o for any value 
of 0o in the range [0, 2tt], we find that the spectrum starts 
as white-noise, i.e. C(k,9o) = A for some constant A, up 
to femin ~ 2ir/L. From that point on, a spurious set of os- 
cillations with period m 4n/L appear superposed on the 
raw spectrum. This illustrated in Fig. |l| in which we plot 
C s (fc, 0). The spurious oscillations merely reflect the 'ring- 
ing' of the Fourier transform due to the sharp edges of the 
window. In the 1-dimensional azimuthal average power spec- 
trum C s (k) — (C s (k,9)) g these spurious oscillations are 
mainly smoothed out, but so too are features in the raw 
spectrum on scales less than A£ m 4n/L (see Fig. l[l]). 
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Figure 2. As in Fig. |l] but for a square multiplied by a cosine 
bell. 



It is well-known (e.g. Press et al 1994) that the spuri- 
ous oscillations in Fourier transforms can be suppressed by 
multiplying the data (sky map) by a window that goes to 
zero smoothly at its edges. In the context of CMB measure- 
ments, Tegmark (1996) calculates the optimal window to 
use in various circumstances to obtain the highest possible 
A£ resolution in the sampled power spectrum. For a square 
field this window is a cosine 'bell' 



W(x) = cos cos 
whose Fourier transform is 



W(k) 



2n cos (k x L/2) cos (k y L/2) 



(vr/L) 



2 JU2 



(10) 



(11) 



At any given point ko in the sampled power spectrum this 
window minimises the variance (second central moment) 
((k — ko) 2 ) of the convolving function |iy(fe)| 2 . Interestingly, 
the more commonly encountered Hann window 



W(x) 



1 + cos ^ 



27ra;\ 



minimises the fourth central moment {(k — ko) 4 } of the con- 
volving function. 

Using the cosine bell window ([lo]), the spurious oscil- 
lations disappear from C s (k, 0) and C(k), but for such a 
small field the features of the raw spectrum are still heavily 
smeared (see Fig. |^) . 

For interferometers we may model the window (primary 
beam) as Gaussian. In Fig. ^| we show the effect on the sam- 
pled power spectrum of a Gaussian window with a w — 0.6 
degrees, which corresponds to a FWHM of 1.3 degrees. Al- 
though for a real interferometer observation, B(k) = 1 only 
where the Fourier domain is sampled (and is zero elsewhere) , 
in order to illustrate the effect of the window we set it equal 
to unity for all k. In Fig. [], we again observe a white-noise 
tail up to fc m in « l/<7m, but from then on the Gaussian 
window does not induce spurious oscillations on the spec- 



Figure 3. (Top) The sampled power spectrum (dashed line) ob- 
served by an experiment with a Gaussian window with a w = 0.6 
degrees, compared to the underlying CMBR spectrum predicted 
by cosmic strings (solid line). (Bottom) The same as above, but 
for the sCDM spectrum. 



trum. The Gaussian window (which occurs naturally for in- 
teferometer observations) is in fact quasi-optimal, and the 
procedure suggested by Tegmark (1996) makes a negligible 
difference to the spectral resolution of the sampled power 
spectrum. However, this window is still too narrow since fea- 
tures in the underlying spectrum on the scale of few x 1 /a m 
are smoothed out. Hence although Gaussian windows do not 
distort smooth spectra, they still erase oscillations in oscil- 
latory spectra if the window is too small. In Fig. |^ we show 
that effect of the Gaussian window on CMBR power spec- 
tra predicted by sCDM and cosmic strings (Magueijo et al 
1995). 

Therefore, for particularly small fields (Ss 2 degrees), 
the sampled power spectrum C s (k) is a bad estimator for the 
underlying CMBR power spectrum C(fc), even after smooth- 
ing out field edge effects. However, in spite of the examples 
given, as the size of the field increases the sampled power 
spectrum rapidly converges to the underlying power spec- 
trum for k > fc m i n . In fact so long as edge effects are dealt 
with, for fields as small as 10 degrees for a square window 
(with cosine bell applied) or 5 degrees FWHM for a Gaus- 
sian window, the sampled power spectrum C"(k) faithfully 
reflects the underlying Doppler peak structure (see for in- 
stance Fig. |E]). For these fields, apart from the white-noise 
tail for k < fc m i n , one then has 

C s (fc) « aC(k)\B{k)\ 2 , 



where the normalisation constant a is given by 
j2 '- r d 2 x „„, n s 



,hk -\w(k)f = 



AW(x)\< 



(12) 



(2tt) 2 ' v " J (2tt) 2 ' v n ~ (2tt) 2 ' 

Equation (Q defines the effective solid angle S7 S of the sam- 
pled field. In the previous plots C(k) was divided by a to 
give the correct normalisation. 
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2.2 Correlations in the sampled power spectrum 

As we have seen, finite fields have the effect of distorting the 
underlying CMBR power spectrum by convolving it with 
the function |VF(fc)| 2 . Aside from smoothing the underlying 
power spectrum, this convolution also has the effect of in- 
ducing correlations between the values of the sampled power 
spectrum C (k) for different values of k. To investigate these 
correlations further, we must first consider the correlations 
induced in the sampled Fourier modes a s (k). 

We expect correlations to exist between neighbouring 
sampled Fourier modes within some correlation length f. 
However, if the window has sharp edges, and the raw C(k) 
spectrum falls-off as 1/k 2 or slower at low k, then in addition 
to nearby correlations there are also 'resonant' long range 
correlations even for relatively large fields. We illustrate 
this point by considering the correlation of sampled modes 
around a given circle in fe-space, i.e. cor[a s (fc, 0), a a (k, 6)\ as 
a function of 6 for some given radius k. This quantity may 
be easily calculated from from (^) or (|^), and is plotted in 
Fig. ^] for k = 50 for various windows. 

For an inteferometer observation with a Gaussian win- 
dow the correlations fall-off at a correlation length of or- 
der f ~ l/cr w around each mode. Furthermore, there 
are no long-range (anti-)correlations, and indeed no anti- 
correlations whatsoever. For single-dish observations with a 
square window, however, apart from correlations between 
neighbouring modes up to f ~ 1/L, there also exist 'reso- 
nant' anti-correlations between modes symmetric about the 
origin. This effect only disappears for k > kt/2, where kt 
defines the angular resolution of the single-dish observation 
as in (|^) . This behaviour may be suppressed by multiplying 
the field by a cosine bell, as discussed above, and the result 
is also shown in Fig. ^. 

The correlations between Fourier modes strongly affect 
power spectra estimation. Let us consider the observed sam- 
pled power spectrum 

f 2w rlf) 

C s (k)= J —\a s (k,6)\ 2 , (13) 

(which in keeping with our notation we write with the s 
downstairs). Although (C s (k)) = C s (k) (as expected), the 
correlation of the Fourier modes induces correlations be- 
tween C B (k) and C a (k') for \k — k'\ < £. These correlations 
are troublesome as they may impart features on any ob- 
served spectrum C s (k) which average out to zero in the en- 
semble average spectrum C 3 (k). This renders C s (k) a bad 
estimator for C(k), and conversely makes C 3 (k) a poor pre- 
diction for what an observers actually sees. For example, as 
we show later, the low k white noise tail C(k) = A (for 
some constant A) corresponds to a single correlated piece of 
the spectrum. As a result a typical observer does not see a 
spectrum C s {k) rs A spectrum with small random fluctua- 
tions around A, but in general each observer sees a spectrum 
of the form C s (k) ~ / A with small random fluctuations 
around fi. Although the ensemble average of jj, equals A, no 
single observer could ever guess it. In general even if the field 
size is large enough so that C{k) is a fair representation of 
the underlying CMBR power spectrum C(k), this does not 
mean that any observed C a (k) would actually let us esti- 
mate C s (k). This, so-called, cosmic covariance problem is 
well known in the context of non-Gaussian theories, where 
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Figure 4. The correlation of sampled Fourier modes 
cor[a s (50, 0), a s (50, 8)] as a function of 9: (top) for a Gaussian 
window with a w = 12 degrees; (bottom) for a square window 
with L = 36 degrees (solid line), and same field but multiplied by 
a cosine bell (dashed line). Correlations fall-off beyond 8 l/a w 
or as 1/L, except for fields with sharp edges, where long distance 
correlations persist. In particular, for a square window there exist 
large anticorrelation between modes symmetric about the origin 
of fc-space. 



it is present even if one assumes all-sky coverage (Magueijo 
1995). 

The best way to deal with cosmic covariance is to do 
away with the correlations. We shall do this by discretizing 
the spectrum estimates in such a way as to obtain a maximal 
finite set of uncorrelated estimates. The separation between 
these uncorrelated estimates will then define the spectral 
resolution Afc(~ A£). 

We could attempt to obtain uncorrelated estimates by 
using the estimator C s (k) as defined in (^) and discretizing 
directly in the space k — \k\. The covariance of C 3 (k) at 
different k values is given by 

f 2lT rlf) 

cov[C s {k),C s {k')] = J —co V [\a s (k,0)\ 2 ,\a s (k',e)\ 2 ],(U) 

which may be simplified significantly if the window W(x) 
used in the experiment is real and even (as are all those 
considered here). In this case its Fourier transform W(k) is 
real, and the real and imaginary parts of the a s (k) modes are 
uncorrelated. For two (possibly correlated) complex Gaus- 
sian random variables z\ and z 2 with uncorrelated real and 
imaginary parts, it can be shown that cov[|zi| 2 , \z2\ 2 ] = 
(ziz^) 2 + (ziz 2 ) 2 , and so 

cov[|a s (fc)| 2 , \a s (k')\ 2 ] = (a s (k)a* s (k')) 2 + {a a (k)a s (k')) 2 . 

Substituting this result into ( p^ ) and using the fact that the 
sky is real, we find 

f 2n Hf) 

cov[C B (k),C s (k')] = 2 ^{a s (k,0)a:(k',e)) 2 . (15) 
Jo 27T 
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The integrand in this expression may be computed from 
(§) or @; an d it can be checked numerically that there is 
always a significant correlation for \k— k | < f, where £ is the 
correlation length discussed above. This can be understood 
from the fact that, roughly speaking, each mode a a (k) is 
correlated with neighbouring modes within a circle of radius 
£ in the fe-plane. Since the esitimator C 3 (k) makes use of 
modes in all directions there are always correlations between 
two estimates C s (k) and C s (fc') for \k — k'\ < £. 

It is clear, however, that the separation Afc between in- 
dependent estimates of the sampled power spectrum should 
increase with k since the number of independent Fourier 
modes per unit k increases. Therefore, instead of using the 
estimator C 3 (k) as defined in (^), and attempting to re- 
move correlations by discretizing directly in k = \k\ space, 
we should first discretize the modes in the two-dimensional 
k space. 



2.3 Uncorrelated meshes and the spectral 
resolution Al 

From equations (^) and (^) we see that the presence of 
the window W(x) causes the underlying (two-dimensional) 
CMBR power spectrum to be convolved with |W(fc)| 2 . To 
obtain independent estimates of the sampled power spec- 
trum, we may divide the fe-plane into a set of discrete cells 
so that for the points at the centres of any two cells (or 
infact any points which have the same relative position in 
each cell) we have 



coi[a 3 (k i ),a* s (k j )] < e, 



(16) 



where fc; and kj are the 'centres' of the ith and jth cell 
respectively, and the constant e< 1. More will be said about 
e and how small it must be in the next section. Although 
we may divide the fe-plane in any way we wish, so long 
as ( |lrj ) is satisfied, we choose to discretize in the form of 
a rectangular mesh. Therefore working outwards along the 
axis corresponding to the smallest dimension of the field, we 
find the set of points for which cor[a s (fci), a s (fci+i)] = e, and 
then repeat the process for the perpendicular axis. These 
point then define our rectangular mesh. Normally the mesh 
will not be very different from a square lattice with cell size 
fco ~ 27r/vn s , where f2 s is the effective solid angle of the 
field as defined in ([l2|). 

From Fig. ^ we see that such a procedure is impossible 
for, say, a square window, because of the long range corre- 
lations it induces. In such a case a the field must first be 
multiplied by the cosine bell (|ll|) before the uncorrelated 
mesh can be constructed. Once the uncorrelated mesh has 
been built we may use a new power spectrum estimator 



C s (fe) 



N(k) 



ifcii=fc 



(17) 



where the sum is over all modes in the mesh for which \ki\ — 
k, and N(k) is the number of such modes. From ( |l7| ) we see 
that the sampled power spectrum is estimated only at a 
finite number of well separated values of k, and that these 
estimates are quasi-uncorrelated. Moreover the estimates are 
distributed as a Xjv(k) distribution. For a square lattice in 
which the modes are separated by fco, the values of k for 
which the sampled power spectrum is estimated are those 



where k — fco ynf + n| for integer numbers m and ri2 . The 
number of modes at a given fc is simply the number distinct 
permutations of integers n\ and ri2 for which \Jn\ + ri2 is 
fixed. For instance, for fc = fco we have iV(fc) = 4, whereas 
for fc = 5fc we find N(k) = 8. It is only for fc > (fco/27r) 
that the number of independent modes between fc and fc + 1 
becomes much larger than 1 (since for large fc the number of 
such modes is « 2-nkjk^). From then on a resolution Al ~ 
Afc = 1 is meaningful and C a (k) should be replaced by 



C s {i) 



N(£) 



(18) 



e<\ki\<e+i 



where N(l) is the number of mesh points for which I < 
\ki\ <£ + l. 

In practice it is the fact that we can only make un- 
correlated estimates with a separation Afc that limits the 
spectral resolution. Spectral leakage also constrains Afc but 
in the context of Doppler peak detection we will see that 
whenever there are enough independent modes to resolve 
the peaks, spectral leakage is no longer a problem. Also in 
that regime the underlying C(fc) spectrum does not change 
by much across each of the cells of the uncorrelated mesh. 
Hence each of the independent estimates C a (k) provided by 
the mesh for fc < (fc 2 /27r) could be placed anywhere in be- 
tween its two neighbours (but of course only at one such 
point at a time). 

Spectral resolution understood in this sense leads to 
the following picture for an experiment with a sky cover- 
age area Q s (for which fc « 2ty/V^s)- The region in the 
spectrum between fc = and fc = fco receives only one in- 
dependent estimate, which bears little resemblance to the 
true underlying spectrum. For fc ^ fco the sampled spec- 
trum and the raw spectrum are roughly proportional, and 
the average separa tion Afc b etween uncorrelated estimates 
goes like Afc m \Jk? + fco fit — fc. For fc ^> fco this means 
Afc « fco/(27rfc), so that the spectral resolution goes as 
1/fc. A meaningful resolution of Al m Afc = 1 can only 
be achieved for fc 3> ko/(2w). For example, if an experi- 
ment has a field of size 10 degrees, then fco ~ 36, and so 
we expect C(k) oc C(fc) for fc <; 60. However, one may not 
achieve independent estimates of the spectrum with sepa- 
ration Al « Afc = 1 before fc 3> 250. Uncorrelated meshes 
provide the maximal spectral resolution Al for a given ex- 
periment. Of course these uncorrealted estimates may be 
grouped into wider bins, which deliberately reduces the spec- 
tral resolution, but also reduces the cosmic/sample variance 
in each of the estimates. 



3 COSMIC VARIANCE IN SMALL FIELDS 

We now compute the cosmic/sample variance in the esti- 
mates provided by uncorrelated meshes both with and with- 
out instrumental noise. Once these are computed, the cos- 
mic/sample variance in estimates obtain by any other bin- 
ning of uncorrelated estimates (as the one proposed in Sec- 
tion is straightforward. 

3.1 Variance of estimators without noise 

Ignoring instrumental noise a centred estimator of the un- 
derlying CMBR power spectrum C(fc) is given by 



© 1996 RAS, MNRAS 000, 



Observability of secondary Doppler peaks 7 



C(k) = 



Cs(k) = 1 
a\B(k)\ 2 N{k)a\B(k)\- 



ifeii=fc 



where a = Sl s /(27r) 2 as defined in equation It is 

straightforward to verify that indeed (C(fc)} = C(k), and 
by performing a calculation similar to that following equa- 
tion ([lij), we find the variance of this estimator to be given 
by 



a\C{k)) = 



2C 2 (k) 



N(k) 

where N(k) is given by 

N(k)=N(k) 2 / J2 {corla^^Uk,)]} 2 



(19) 



(20) 



The quantity Nk acts as an effective number of independent 
modes contributing to the estimate C(k). If all the N(k) 
modes were perfectly correlated or anticorrelated, then they 
would only count as one, and in that case N(k) = 1. If, on 
the other hand, all the N(k) modes are very nearly uncor- 
rected, then N(k) ~ N(k). Following a similar calculation 
to the above the correlations between two estimates can is 
found to be 

coT[C{k),C(k')} = Arn , 2 Arn ,, Y" {cor[o s (fei),a:(fej)]} 2 . 



N(k)N(k>) 



If an 'uncorrelated' mesh is built such that for any two mesh 
points 

cor[a s (fei), a*(fej)] < e, 



as discussed in section 2.3, then residual correlations of order 



^ ' I y N(k)N(k') will persist between the estimates C(fc) and 

C(k'). In this case we also find N{k) « N(k)(l - e 2 ). 

These results allow a more concrete definition of the 
uncorrelated mesh. The point of uncorrelated meshes is to 
do away with correlations among estimates without throw- 
ing away information. However reducing the cosmic/sample 
variance in an estimate and reducing the correlations be- 
tween estimates are contradictory requirements. In other 
words, by lowering the value of e arbitrarily, one reduces 
correlations between the C(k) estimates, but their individ- 
ual variances increase. 

The best compromise is obtained by considering the 
effective number of independent modes within a mesh cell 
centred on ki, which is given by 



iV ce ll(fei) = 1/ 



d 2 k 

A-ell 



{cor[a 3 (fci), a*(k)]} 2 , 



where ^4 cc ii is the area of the mesh cell. Clearly, N cc n(k) is 
always greater than unity. One can therefore, on average, 
avoid the loss of non-redundant information by computing 
the average density of independent modes around a given 
mode at k — ki 

p(ki) = 1/ J d 2 k{coT[a s (ki),a* s (k)]} 2 , 

and defining the mesh size as ko = 1/ »/p. In this way power 
spectrum estimators derived from large regions of the k- 



plane have the same variance whether one uses the mesh or 
the continuum of modes in its calculation. 

In well behaved regions where C s (k) oc C(k) we find 
that that fco = V2~7r/a w for a Gaussian window, and (of 
course) fco = 2n/L encodes all the information in a square 
field (multiplied by a cosine bell). This choice of fco corre- 
sponds to e ~ 0.2. We may may check that even such a large 
value of e is small enough to make all correlations negligi- 
ble. Since in the estimator C(k) no two neighbouring cells are 
ever used for the same fc, one may safely set N(k) ~ N(k) 
when computing cr 2 (C(fc)). Furthermore the residual corre- 
lations between C(fc) and C(fc') involves at most two neigh- 
bouring modes, and so are of the order 0.04/iV(fc), which is 
typically smaller than 1 per cent. Hence one need not throw 
away any information in order to keep residual correlations 
low, as uncorrelated meshes spread the modes among esti- 
mates in the best possible way. 

For fc ^> (fc ( 2 /27r) individual G coefficients may be es- 
timated by 



Cs(£) = 1 
a\B(£)\ 2 a\B(£)\ 2 N(£) 



y M fc i)l 2 



K\ki\<l+1 



For a square window of side L the variance in this estimator 
is given by 



2 [C(£)] 



2C 2 (£) _ C 2 (£) 4tt 
N(£) 



£ L 2 ' 



(21) 



since N(£) w N(£) w 27r^/fcg » L 2 l/2n. Equation @ has a 
simple interpretation. The cosmic variance of C(£) is simply 
2C 2 (£)/{2£ + 1) w C 2 {£)/£ for large £, and 4tt/L 2 is the 
fraction of the sky covered by the experiment. However, this 
formula is only a good approximation for £ w fc 3> (k%/2ir), 
and significant corrections are necessary for smaller values 
of I. 

Finally, we note that by increasing the size of the ob- 
served field one does not decrease the relative error in the 
estimate provided by each of the modes a a (k). The relative 
error in each of the modes is fixed by its Gaussian nature, 
and its value <r 2 (|a s (fc)| 2 ) = 2C s2 (fc) merely reflects the van- 
ishing kurtosis of a Gaussian distribution. The positive effect 
of increasing the size of the field is that it populates Fourier 
space more densely with uncorrelated modes. This improves 
the spectral resolution and also reduces the cosmic variance 
in the estimates based on a larger number of independent 
modes. 



3.2 Observations of multiple fields 

So far we have assumed that observations are of a single fi- 
nite field. It may sometimes happen, however, that the total 
set of observations consists of rif > 1 fields each observed 
with a window Wi(x) for i = 1 to n/. If a s (k) and b a (k) are 
the sampled modes derived from observations of the ith and 
jth field respectively, then it is straightforward to show that 
their covariance matrix is given by 



(a s (k)b* s (k')) = 



d 2 k'' 
(2,r) : 



C(k")Wi(k - k")Wj{k' - fc").(22) 



If the jth window can be obtained from the ith window by 
a translation with vector R then 
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Wj{k) = Wi{k)exp(iR- k), 

and the only chance for a correlation is therefore if k — 
k' . However, whenever the approximation C s (k) oc C(k) is 
valid, we find that 



cor[a s (fc),6*(fc)] 



J^\W(k)\ 2 exp(ik-R) 



I 



d 2 k 



\W(k)\ 



J d 2 x W(x)W(x 



R) 



J d 2 x\W(x)\ 2 

which is of the order of the percentage of area overlap be- 
tween the windows. Therefore, if the fields are well-separated 
then correlations between the sampled modes from each field 
are very small. In such a case one can build an uncorrelated 
mesh for each field, and superpose the meshes giving an 
extra index i = to each of the mesh points. The 

uncorrelated mesh estimates in well behaved regions [where 
C(k) oc C(k)] will then have the same properties as the 
single field estimates but with N(k) replaced by rifN(k). 
Therefore the relative error in the uncorrelated estimates of 
the sampled power spectrum can be reduced by observing 
several well-separated fields, or (for single dish observations) 
by dividing observations of a single large field into several 
smaller (but non-overlapping) fields. 



4 ESTIMATING THE CMBR POWER 

SPECTRUM FROM REAL OBSERVATIONS 

So far we have only considered the effects of the window and 
the beam of an experiment on the sampled power spectrum. 
Real observations, however, have the added difficulties of 
instrumental noise, the discretization of the sky temperature 
distribution into pixels (for single-dish experiments) , and the 
existence of foreground emission. 

4.1 The effects of noise and pixelisation 

We first consider the effects of noise and pixelisation on 
sampled power spectrum estimates drawn from uncorrelated 
meshes. In our analysis we neglect residual correlations as 
well as any terms of order e 2 . 

Let us first consider single-dish observations, in which 
the CMBR sky map of the field is discretized into M pixels 
at positions Xi and with noise Ni for i = 1 to M. For simplic- 
ity we shall assume that the Ni are independent Gaussian 
random variables with (Ni) = and 

This assumption is a reasonable approximation to reality. 
To simplify our discussion still further we also assume that 
the variance of the noise in each pixel is constant so that 
<y% = Opj x . We can eliminate the discreteness by instead con- 
sidering the noise as a continuous Gaussian random field 
N(x) which is added to the sky temperature distribution, 
and is characterised by the two-point correlation function 



{N(x)N(x')) « 6(x - 



7 pix^pix 



(23) 



where f2 p i x is the area of a pixel. Assuming all-sky cover- 
age these two quantities are often combined into w^ 1 — 
(jp ix r2pi x , where w may be considered as a weight per unit 



solid angle (Knox 1995). If ttot is the total observational 
time available, then the time spent observing any given 
pixel is tpix = itotfipi x /(47r). The noise variance per pixel 
°pix — s 2 A P ix, where s is the sensitivity of the detector. 
Therefore, for a detector with fixed sensitivity s and for a 
given total observation time ttot, by changing f7 p i x the quan- 
tity u) _1 = 0"p ix f2 p ix = Aits 2 /ttot remains constant. Hence 
io _1 is important as a qualifier for noise in all-sky maps ob- 
tained using different scanning strategies. If we consider the 
most general case where we make rif maps each of size Q s , 
then the time spent on each pixel is t p i x = ttot^pix/(^f^s), 
and the the quantity which now remains constant is 

W ■ ' (24) 



( 



ttot p>*--t— \ n{ Q g 

The noise defined by (^3|) adds an extra term in the 
covariance matrix of the sampled modes (a s (k), a s *(k')) of 
the form 

{a s {k)a a (k )} = v / d x\W(x)\ e . (25) 

Since an actual sky map is divided into pixels of size d (say), 
we would estimate the sampled Fourier modes a s (k) using a 
Fast Fourier Transform (FFT). We would therefore obtain 
estimates of the Fourier modes between k — and k — 2iv/d 
on a square grid with spacing 2n/L, where L is the size of the 
field, i.e. one estimate per uncorrelated mesh cell. Therefore, 
neglecting terms of order e 2 in (53), we find that the noise 
on the uncorrelated mesh points is approximately 

(a„ (ki)a s *(fcj)> = Jyacrpixfipix, (26) 

where a = Q. s /(2ir) 2 as defined earlier. 

For interferometers, on the other hand, instrumental 
noise is added in Fourier space directly, and may be quanti- 
fied by means of the noise per mesh cell a%, such that 



(a s (ki)a s *(kj)) = S tJ a N 



(27) 



We shall assume that the interferometer possesses sufficient 
antennas appropriately positioned to provide roughly uni- 
form coverage of the Fourier domain in the fc-range of inter- 
est. A method for positioning a given number of antennas 
in such a way is given by Keto (1996). It may, of course, 
be possible to tailor the coverage of the Fourier plane so 
the density of measured visibilities is greatest in the regions 
of the (one-dimensional) fc-range which are of most interest 
(so that the interferometer may be thought of as a matched 
filter to the expected CMBR power spectrum). This does, 
however, rather preempt the results of the observation, and 
furthermore, in addition to power spectrum estimation, in- 
terferometer experiments are usually designed to make maps 
of the CMBR, for which uniform coverage is desirable. 

If the density of measured visibilities (i.e. the number 
of visibilities per unit area in the fe-plane) is p v is, then the 
number of visibilities per mesh cell is simply 



1.2- 27r 
^vis — Pvis"-0 ^ Pvis p. 



where ko is the cell size. The measured visibilities in each cell 
(say at fc = fci) may then be used to estimate |a s (fci)| 2 by, 
for example, a maximum likelihood analysis (see Hobson, 
Magueijo & Kaiser, in preparation), or by some weighted 
average in the limit of low noise. In any case, a lower limit 
to the noise per mesh cell is given by 
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2 s A n 

a N = 



(28) 



where i v j s is the time spent measuring each visibility. The 
presence of Q 2 (the area of the observed field) in the first 
equality results from defining the Gaussian window of the 
experiment to be always equal to unity at its peak. By chang- 
ing the spacing of the inteferometer antennas we may fix the 
limits fe m in and fc max of the fc-range to which the experiment 
is sensitive. In targetting differences between the cosmolog- 
ical theories, or determining the existence (or otherwise) of 
Doppler peaks, these limits are usually fixed. In any case, 
the the area of the fc-plane sampled is A = 7r(fcn lax — fc^m)- 
Now the time spent observing a particular field is given by 
t{ = Apvistvis, but the total observing time t to t is divided 
among nf fields then the time per field is simply tf — t tot /ri{. 
Therefore for comparing inteferometer observations with the 
same detector sensitivity s 2 , total observing time ttot and fc- 
plane coverage A, but in which the size of the field f2 s and 
the number of such fields ri{ can vary, then the quantity to 
keep constant is 



_! As' 
w = 

ttot 



We may then describe the noise per mesh cell by 



(29) 



{a s (ki)a s *(kj)) = 5 ir 



Noise always introduces a bias in the estimator C(k). 
Hence in the presence of noise one should use instead the 
centred estimator for single-dish observations 



C(k) = 



n f N(k)a\B(k)\ 2 ^ 

|fcj | — k 



2 
°pix i2 pix 

WW 



(30) 



and for interferometers 
f 



C{k) = 



n{N(k)a ^-^ 1 v " / a 

\ki\=k 



(31) 



For k S> (/cq/27t) these estimators should be rewritten in the 
usual way. The variances in these estimators are increased 
by the presence of noise. For single-dish experiments 

2 



C 2 (k) 



2 f <x»plx V 

n f N(k) [ \B{k)\ 2 C(k) J 



1 + 



w Q 3 nf 



n f N{k) \ 4n\B{k)\ 2 C(k) 
whereas for interferometers 



(32) 



v 2 (C(k)) 
C 2 (k) 



n { N{k) 
2 

n t N{k) 



1 + 



1 + 



aC(k) 



C(k) 



(33) 



4.2 The effects of smooth foreground emission 

Depending on the frequency at which the observations 
are made, the existence of various types of foreground 
emission can severely hamper the measurement of CMBR 
anisotropics. A discussion of these foreground components, 



and the regions of frequency /multipole space in which each 
dominates, is given by Tegmark & Efstathiou (1996) (TE96). 
The main components of this foreground emission are extra- 
galactic radio point sources, and continuous Galactic dust, 
synchrotron and free- free emission. 

The general method for dealing with continuous Galac- 
tic foreground emission is to make observations at several 
different frequencies, and use this spectral information (to- 
gether with their predicted power spectra in the Wiener fil- 
tering approach as in TE96) to perform a subtraction of 
these components. This process results in estimates of tem- 
perature maps or Fourier modes due to the CMBR alone, 
together with some errors on these estimates. 

In general the power spectrum of the these separation 
errors is not strictly constant across the fc-plane, but may 
often be approximated as such. We may then model the re- 
construction errors by the presence of a generalised noise 
field, and include the effects of foreground subtraction sim- 
ply as a factor which multiplies the noise per mesh cell. 
Therefore, in the expression (for single-dish observations) 
for ui -1 (p4J) , the quantities <T p j x and f2 p j x refer to the pixel 
noise and pixel solid-angle in the final CMBR temperature 
map deduced from observations at several frequencies using 
some separation algorithm. Similarly, crjv in p9[ ) refers to 
the error on the deduced CMBR Fourier modes. 

Clearly the size of the errors in the CMBR map or 
power spectrum will depend on the value of ui _1 for each 
frequency channel, on the number of such channels, and on 
the separation algorithm used. TE96 show that for a 'pixel- 
by-pixel' subtraction algorithms these errors are typically 
several times larger than the average error on an individual 
frequency channel due to instrumental noise alone. Using 
a vector Wiener filtering algorithm, however, these errors 
can be reduced by up to a factor of 10, depending on the 
number of frequency channels and their range in frequency. 
Nevertheless, as the authors themselves suggest, the Wiener 
filtering method is rather optimistic in that it assumes that 
both the frequency dependence and the power spectra of the 
various foreground components are known reasonably well, 
which is certainly not the case. A reasonable estimate might 
be that the errors associated with the separation process are 
of a similar magnitude to the average errors on an individual 
frequency channel due to instrumental noise alone. 

Finally, we note that since we are using an uncorrelated 
mesh in the estimation of the CMBR power spectrum, the 
additional errors on the power spectrum estimates arising 
from any foreground separation process will still be uncor- 
related. 



4.3 The effects of point sources 

As mentioned above the existence of foreground point 
sources can cause problems for estimating the CMBR power 
spectrum. In general point sources cannot be removed from 
spectral information alone, but require the identification of 
the sources by higher-resolution observations. The general 
scheme is to survey the same region of sky as that observed 
by the CMBR experiment at a frequency close to that of 
the CMBR observations in order to indentify all the point 
sources down to some flux limit, which are then subtracted 
from the CMBR data. The required flux limit is usually 
such that the confusion noise from unsubtracted sources is 



© 1996 RAS, MNRAS 000, 



10 M.P. Hobson, Joao Magueijo 



roughly equal to the instrumental noise. Once this source- 
subtraction has been performed the residual emission from 
unsubtracted sources can be modelled as a an additional 
smooth foreground component, and treated in the same way 
as the Galactic emission, as discussed above. 

We note here that although it is generally believed that 
point source contamination becomes less important as the 
observing frequency increases above about 100 GHz, there 
is no direct evidence for this. Moreover, even the popula- 
tion of radio point sources at frequencies above about 10 
GHz is rather uncertain, and it may be inadvisable to rely 
on low frequency surveys such as the 1.5GHz VLA FIRST 
survey (Becker et al. 1995) to subtract point sources from 
CMBR maps made at much higher frequencies. In order that 
one can be confident in the final CMBR map/power spec- 
trum, it is therefore necessary to make higher-resolution ob- 
servations at frequencies close to those of the CMB experi- 
ment. Therefore, if the CMB observations are made over a 
large frequency range, higher-resolution observations at sev- 
eral frequencies in this range may be required. Clearly, this 
can severely limit the maximum possible sky coverage, and 
this should be borne in mind when interpreting the figures 
presented in the next section. The detailed implications of 
source subtraction for total sky coverage will be discussed 
fully in a forthcoming paper. 







/ 
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Figure 5. sCDM (line) as seen by a square window with L = 10 
degrees after and multiplication by a cosine bell (dash). The 
points indicate the corresponding power spectrum estimates for 
a particular realisation in each of the bins indicated by the hori- 
zontal error bars. The vertical error bars indicate the theoretical 
cosmic/sample in the absence of noise. We have set B(k) = 1 for 
clarity. 



5 OBSERVING SECONDARY PEAKS OF 
STANDARD CDM 

If one is interested only in determining whether or not sec- 
ondary oscillations exist in the CMBR power spectrum, then 
the uncorrelated mesh estimates should be combined into 
broader bins, centred on the anticipated positions of the 
peaks and troughs in the power spectrum under consider- 
ation. Six broad bins adjusted for testing standard CDM 
(sCDM) are shown in Fig. ^. The figure shows the ensemble 
average sampled power spectrum for an experiment with a 
square window with L = 10 degrees (dashed line), multi- 
plied by a cosine bell, as compared to the underlying power 
spectrum (solid line). The points indicate the average power 
in each of the bins denoted by the the horizontal errorbars. 
These points are calculated by averaging the values of C(k) 
at each measured A;-value in each bin. The vertical error- 
bars are due to cosmic/sample variance on these estimates 
(assuming no noise). They were computed by adding the 
variances of the measured C(k) estimate values in each bin, 
and dividing by the square of the number of such estimates 
in the usual manner. 

Once the bins have been chosen a statistic for the pres- 
ence of oscillations in the CMBR power spectrum can be 
derived by inferring the convexity of the power spectrum at 
each bin position (apart from the first and last bins). If the 
estimated power in the ith bin in C\ (where here i = 1 to 
6), then we define (for i = 1 to 4) the quantities 



Ci 



d-i + Ci-i 



a. 



(34) 



These 'convexities' are all negative if there are no secondary 
peaks, but alternate in sign for sCDM. The variance in the 
ith convexity is simply 



a 2 (d) = a 2 (d) + 



a 2 (Ci-i) + a 2 (Ci+i) 



(35) 



Therefore a convexity in bins 2 or 4 which is larger than zero 
by an amount n x a(d) corresponds to a n-sigma detection 
of secondary oscillations. Hence we can define an oscillation 
detection function as 

v _ \(&)\ 



c{Ci 



(36) 



for the relevant i — 2 and i = 4. The function E, tells us 
to within how many sigmas we can claim a detection of sec- 
ondary peaks. More generally we could define in a similar 
way a set of oscillation statistics centered at a generic point 
£ — £i in the power spectrum and making use of bins with 
width A£. These statistics would then measure the oscilla- 
tory character of the spectrum as a function of the point 
li and the scale A£. In what follows, however, we shall al- 
ways target particular theories, and so we will not use this 
broader class of statistics. 

In Table [l] we display the results of applying the convex- 
ity statistic to sCDM under the experimental situation used 
in Fig. ^| and we see that a very significant detection can 
be obtained for this relatively small field (L = 10 degrees) 
in the idealized case of no noise (w" 1 = 0) and infinite res- 
olution (ai, — 0). In Fig. ^ we plot the dependence of the 
detection functions E2 and E4 on sky coverage for the same 
idealized case. From the figure it is clear that a total sky 
coverage of greater than about (5 deg) 2 is necessary for a 
detection of secondary peaks. This supports our earlier claim 
that the convolution of the power spectrum is never a prac- 
tical problem in this context (once the field edges-if they 
exist- have been smoothed by a multiplication with a cosine 
bell). The deterioration of the detection function for fields 
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140 


300 


5.21 


0.78 






330 


420 


2.69 


0.54 


1.67 


3.75 


430 


570 


3.53 


0.51 


-0.73 


2.25 


570 


690 


2.89 


0.46 


0.64 


2.27 


700 


820 


3.54 


0.50 


-1.26 


4.61 


900 


1050 


1.67 


0.30 







Table 1. Bins (£1,(2) used if Fig. g, the average power in bin d 
and its variance, and convexities d together with the number of 
sigmas within which one could be certain of observing their 
correct sign. 
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Figure 6. £2 (line) and £4 (dash) functions for interferometers 
(top) and single-dish experiments (bottom). We have assumed no 
noise (w -1 =0) and infinite resolution (05 = 0). For the top plot 
9 is the FWHM of the interferometer primary beam in degrees 
(FWHM « 2.35ct„). For the bottom plot 6 is the side L of a 
square field in degrees. 

smaller than this is, however, considerably worse than sim- 
ple extrapolation of our curves in Fig. |i| imply. From Fig. ^ 
we also see that the first dip in the sCDM power spectrum is 
more easily detected than than the second one, a situation 
only exarcebated by finite resolution and the prescence of 
instrumental noise. Therefore, in the rest of this section, we 
shall confine ourselves to considering the detection function 
E2, which from now on we refer to simply as E. 

We now undertake a general analysis of the E function 
for single-dish and interferometer experiments separately, 
taking into account the effects of instrumental noise and 
finite resolution. 

5.1 Single-dish experiments 

For single-dish observations the detection function E = 
E(w _1 , 8b, L), i.e. it is a function of the noise level u> -1 
(which itself depends on the total observation time and 




Figure 7. Contours £ = 1, 2, . . . , 10 for an all-sky experiment. 
The noise ui _1 varies from (5/iftT) 2 (deg) 2 to (125/i.fQ 2 (deg) 2 . The 
no-noise axis (not plotted) is the 76.6 contour which asymptoti- 
cally absorbs all other contours. 

detector sensitivity), the beamsize characterised by 0b ~ 
2.35<7(,, and the size of the observed field L. We now con- 
sider two sections through this (three-dimensional) function 
for fixed L. 

The first section is the all-sky limit, which may be ob- 
tained in our formalism, not by considering an infinite patch 
of sky, but by setting L 2 = 4ir. This is equivalent to setting 
L ~ 202 degrees. The all-sky section, plotted in Fig. 0, shows 
first of all that the main consideration for detecting the sec- 
ondary peaks is not the beam size, but the noise level. In 
the absence of noise = 0) one could in theory decon- 

volve a beam of any size without adding any extra uncer- 
tainty. It can be checked that in the all-sky case the w" 1 — 
axis is the contour E ~ 77 of the detection function. This 
is also the maximum possible detection, imposed by cos- 
mic variance, for an an all-sky single-dish experiment (with 
any beam). As soon as noise is taken into account, however, 
this upper limit is greatly reduced, and the beamsize be- 
comes crucial. Even for the low (but realistic) level of noise 
w^ 1 = (15/iK) 2 (deg) 2 a modest 1-sigma detection requires 
6b = 0.6°, whereas a 3-sigma detection requires 0b = 0.5°. 
This openly contradicts the naive argument favouring any 
beam with 0b < 0.85°. Although such beams do not start to 
cut off before I = 1000, the noise is the real limiting factor. If 
the level of noise is pushed up to to -1 = (125/iK ) 2 (deg) 2 an 
all sky experiment would never provide more than a 2-sigma 
detection, and even this would require 0b = 0.1°. 

The all-sky diagram is misleading, however, as for some 
levels of noise and beam sizes it may be advantageous to 
reduce the sky coverage. In Fig. |^ we plot the detection 
function for L — 20 degrees. We see that all combinations 
of noise and beamsize outside the E = 6 contour for an 
all-sky experiment profit from reducing the sky coverage. 
We call this the noise-dominated region, and its comple- 
mentary region signal-dominated. The effects of noise are 
always reduced by concentrating all the observation time on 
a smaller area, since then we allow the same coherent sig- 
nal to compete with the incoherent noise. It therefore makes 
sense that, if one is in the noise-dominated region, the detec- 
tion function increases by reducing the area of sky observed. 
In the signal-dominated regions, on the other hand, the de- 
tection function is controlled by the cosmic/sample variance. 
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Figure 8. Same as Fig. [?j, but for a square field with L = 20 
degrees. 

Different observers in the Universe looking at different sky 
patches see different signals scattered about the theoretical 
average with a variance equal to the cosmic/sample vari- 
ance. Although we cannot increase the number of skies ob- 
served, we can increase the area of sky sampled, thereby 
decreasing the cosmic/sample variance. For a given sample 
size the cosmic/sample variance limit is achieved at the no- 
noise axis and is roughly E « 77 L /(47r). By increasing L 
the cosmic/sample variance decreases and so the detection 
function increases in signal-dominated regions. 

The contrasting behaviour of the noise- and signal- 
dominated regions can be explained by the following semi- 
quantitative argument. From Eqn. (p3) we have roughly 

« 7 (37) 

where d and Bi are the average power spectrum and beam 
value for the bins used in estimating the convexity E. For a 
low signal-to-noise ratio, the noise term in the denominator 
dominates and so E 2 oc 1 /f2 s . On the other hand, for a high 
signal-to-noise ratio, the noise term is much smaller than 1, 
and so E 2 oc 1/Sl* - /3fT 2 , where /3 < 1. 

The question naturally arises as to what is the ideal 
scanning strategy for a given experiment, and assuming that 
ideal strategy how significant is the detection of the sCDM 
secondary Doppler peaks. The answer to this question de- 
pends on the values of noise level w~ 1 , which itself depends 
on current detector technology and how effectively fore- 
ground emission can be subtracted. If we assume a low noise 
value (such as, for definiteness, w -1 = (25/iLT) 2 (deg) 2 ), then 
the contours of E in the beam-sky coverage plane are as 
shown in Fig. || For any beamsize there is a maximum sky 
coverage beyond which the detection is not improved. If any- 
thing the level of the detection decreases, but typically not 
by much. The ideal scanning strategy is then defined by a 
line Li(8t) which intersects the contours of E at the lowest 
L-value at which a plateau has been achieved in the detec- 
tion function. The significance of the detection obtained for 
an ideally scanned experiment depends on the beam size. 
For example, if 9b = 0.6°, the ideal coverage is a patch of 
Li (0.6°) = 5 degrees, which results in a 3-sigma detection. 
If 9i — 0.5°, on the other hand, an 8-sigma detection can 
be obtained with Li — 35 degrees. The detection provided 




L (degrees) 



Figure 9. Low noise (ui _1 = (25/ii£~ ) 2 (deg) 2 ) contours of the 
£ function. We allow the sky coverage L to vary between 5 de- 
grees and all sky (L = 202 degrees), and consider beamsizes with 
FWHM between and 1 degrees. 
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Figure 10. Very low noise (w^ 1 = (\5fiK) 2 (deg) 2 ) contours of 
the £ function. 



by an optimally scanned experiment increases at first very 
quickly as the beam is reduced below 9b = 0.6° (from 3- 
sigma at 9 = 0.6° to 33-sigma at 8b = 0.2°). By reducing 9b 
from 0.2° to zero, however, the detection is only increased by 
2-sigma (from 33 to 35). For this level of noise the maximal 
detection is 35 sigma and is achieved with 9b < 3' and all- 
sky coverage. For low noise levels all-sky coverage is never 
harmful, but it is the beamsize that determines how good 
a detection can be achieved, and how much sky coverage is 
actually required for an optimum level of detection. 

For noise levels of the order w' 1 = (25^L") 2 (°) 2 the 
overall picture is always as in Fig. ^. In particular, there is 
always a top contour (like the E = 35 contour in Fig. ^ re- 
ferring to the maximal detection allowed by the given noise 
level. The maximum E is always achieved with infinite reso- 
lution, but one falls short of this maximum by only a couple 
of sigmas if 9b ~ 0.1°. If the noise is much smaller than this, 
however, the summit of E is beyond L = 202° , as in Fig. |To[ 
For w~ x = (15 /J-K) 2 (deg) 2 , for instance, all-sky coverage be- 
comes ideal for any 9 b < 0.3°. 

If, on the other hand, the noise is much larger than 
= (25(iK) 2 (dee) 2 then the E contours are qualitatively 
different from Fig. U. For, say w~ - 1 = (60 f iK) 2 (deg) 2 , the 
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Figure 11. High noise (w~ 1 = (GO^iK) 2 (deg) 2 ) contours of the 
£ function. 

contours of E are shown in Fig. |ll]. The beamsize is now a 
crucial factor. A beamsize of 9b = 0.5° would provide a 3- 
sigma detection (with Li — 10°), but reducing the beamsize 
to about 6b = 0.4° improves the detection to 6-sigma (with 
Li = 20°). Figure Fig. jTl] also shows that, for high noise 
levels, forcing all-sky coverage dramatically decreases the 
detection. 

We note finally that in interpreting the figures presented 
in this section, we must remember that in order to achieve 
low values of the necessity for point source subtrac- 

tion may severely limit the maximum possible sky coverage. 
In turn this can greatly reduce the significance of possible 
detections. 



5.2 Interferometers 

For inteferometer observations the detection function E = 
E(tu -1 , 6 W , rif), so it is a function of the noise level the 
size of a single field (which is characterised by the FWHM of 
the Gaussian primary beam 6„ « 2.35a w ), and the number 
of such fields nf. Alternatively, noting that the solid-angle of 
a single field Q, B = ira^, we can instead write the detection 
function as E = E(w~ ,6 w ,ilt), where fi t = rifQ B is the 
solid angle of the total area of sky observed. 

For ground-based interferometers it is expected that at- 
mospheric effects will limit the shortest possible antenna 
spacing, which in turn will restrict the largest allowable 
FWHM of the primary beam to be ~ 5 degrees. Since this is 
also close to the limit where convolution of the power spec- 
trum becomes important for sCDM, in this section we shall 
restrict our attention to this case. 

In Fig. ^we have plotted the contours of E for 8„ = 5°. 
As in the single-dish case, we we find that for any given noise 
level there is a sky coverage beyond which the detection 
function saturates and then starts to decrease. As the noise 
level decreases, the the saturation sky coverage increases 
along with the significance of the detection achieved. 

In the high-noise regime w -1 > (50fiK) 2 ra,d~ 6 the ideal 
scanning strategy is to observe only a small sky area (a 
few fields at most), which results in only a 1- to 2-sigma 
detection. Increasing the sky coverage in the high-noise 
regime seriously decreases the level of the detection. We note 
that a 1-sigma dectection requires w _1 < (WOfiK ) 2 rad -6 , 




sqrt(sky coverage) (degrees) 



Figure 12. Contours of the detection function £ for an inter- 
ferometer with a primary beam with 8 W = 5°, for varying noise 
levels ui _1 and sky coverage. 

and a 2 sigma detection may be achieved with w^ 1 = 
(50fiK) 2 rad~ 6 . For w' 1 < (50/iJf) 2 rad -6 we enter the sig- 
nal dominated region. By reducing w -1 both the signifi- 
cance of the detection and the ideal sky coverage (no. of 
fields) increase very quickly. We note that observing just a 
single field provides only a 2-sigma detection dictated by cos- 
mic/sample variance. By increasing the number of observed 
fields up to a saturation value, however, increases the level of 
the detection significantly. In the signal-dominated regime, 
increasing the number of fields beyond the ideal value de- 
creases the detection only very slightly. For low noise levels 
w (W/iK) 2 rad~ 6 we see that an 9-sigma detection is 
possible for an ideal total sky coverage of fi t ~ (45 deg) 2 . 
We also note that in the low-noise case, detections of sev- 
eral sigma are possible even for relatively small sky-coverage. 
This is an important consideration if source subtraction is 
taken into account, since (as mentioned earlier) this severely 
limits the total possible sky coverage. 



6 DISTINGUISHING BETWEEN COSMIC 
STRINGS AND INFLATION 

It is shown in Albrecht et al (1995) and Magueijo et al (1995) 
that the power spectrum of cosmic strings does not possess 
secondary Doppler peaks. It instead has a single peak at 
multipoles around £ = 400 — 600, although the height and 
precise shape of the peak are still subject to a considerable 
theoretical uncertainty. 

If future observations of the CMBR show that the main 
peak in the true CMBR power spectrum is at lower multi- 
poles, such as I = 200 — 250 (as predicted by sCDM), then 
cosmic strings (in their present form) can be rejected as a 
possible theory for structure formation in the Universe. If, 
however, the main peak of the true CMBR power spectrum 
is shifted towards higher multipoles, then we must rely on 
the presence or absense of secondary peaks in order to distin- 
guish between cosmic strings and CDM/inflation scenarios. 

In this section we consider the case of maximal confu- 
sion by comparing a cosmic strings model and a CDM model 
for which the main peak in the power spectrum has the same 
position and shape (but the latter exhibits secondary peaks). 
For definiteness we have chosen a CDM theory with a flat 
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Figure 13. The angular power spectrum of stCDM (solid line) 
and one possible cosmic string scenario (dotted line). The dashed 
line is the stCDM power spectrum C s (k) as sampled by an in- 
terferometer with a primary of 2 degrees FWHM. The points in- 
dicate the average power in each bin for stCDM. The horizontal 
errorbars denote the width of the bins, and the vertical show the 
sample variance of the power estimtates for the same interferom- 
eter, assuming no noise. 



primordial spectrum, 0, — 0.3, h — 0.6, and Qth 2 = 



0.02. 



We shall call this theory stCDM, the CDM competitor of 
cosmic strings. We have normalised the power spectra so 
that the height and shape of the main peak are similar in 
both cases. In principle we could now study the convexities 
associated with bins adjusted to the stCDM secondary os- 
cillations, study the same quantities drawn from a cosmic 
strings spectrum, and then calculate the difference between 
the convexities average values in the two theories in units of 
the average variance. However, this procedure would define 
a 'distance' between the two theories which depends on the 
shape of the strings peak. We therefore simply study the 
first dip detection function of stCDM, and then take this 
detection function as a cosmic string rejection function. 

In Fig. [l3] we show the angular power spectrum of 
stCDM (solid line) and a possible power spectrum for cosmic 
strings (dotted line). We also plot the stCDM power spec- 
trum C s (k) (divided by a) as sampled by an interferometer 
with a primary beam of 2 degrees FWHM (dashed line) . We 
see that, since the secondary peaks for stCDM are now much 
more separated in £ space (as compared to sCDM), the field 
size required for convolution of the power spectrum not to 
be a problem is much smaller. It can be checked that inter- 
ferometers with a primary beam of FWHM <; 1.5 degrees, 
and square windows (multiplied by a cosine bell) with L J; 4 
degrees are perfectly acceptable. We therefore assume these 
values as lower bounds on the field size for the purposes of 
cosmic string detection. For fields smaller than this decon- 
volution of the power spectrum would be necessary causing 
a deterioration of the detection function not taken into ac- 
count in our calculations. However, we will see that the field 




Figure 14. Low noise w 1 = (25[iK) 2 (deg) 2 contours of the 
detection function £ for stCDM. 




Figure 15. High noise w 1 = (QOfiK) 2 (deg) 2 contours of the 
detection function E for stCDM. 



must always be larger than this for any reasonable detection 
to be achieved, even baring deconvolution effects. 

In Fig. hj we also show the bins in £ chosen for the 
study of the convexity associated with the first dip in the 
stCDM power spectrum. The vertical errorbars are the sam- 
ple variance errorbars associated with these bins for an in- 
terferometer with a primary beam of 2 degrees FWHM. We 
then repeat the same exercise as in the previous section con- 
cerning the detection function of the first dip of stCDM. In 
Figs. |lj, |l^, and |l6| we redraw the experimental parameter 
space covered in Figs. |, 0, and|l| but as applied to the first 
dip stCDM detection function. The only difference is that 
now we allow L to start at 2 degrees. Similar, for interfer- 
ometers, we consider a primary beam of 2 degrees FWHM. 
Again we make the point that if the effects of source subtrac- 
tion are taken into account the total possible sky coverage 
may be severely limited, and this should be borne in mind 
when considering the figures presented in this section. 

Overall we see that in signal dominated regions the de- 
tection is much better for stCDM than for sCDM. This is be- 
cause features at higher £ have a smaller cosmic / sample vari- 
ance (which is proportional to l/l). It can be checked that 
the cosmic/sample variance limit, obtained with a single- 
dish experiment with no noise, is now E ~ 197L 2 /(47r) (as 
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2° fields, in order to obtain a detection between 1- and 2- 
sigma. For w" 1 < (l50fiK ) 2 rad -6 we enter the signal domi- 
nated region. Following the ideal scanning line with decreas- 
ing ui _1 , both the sky coverage and the level of the detection 
start increasing, at first first slowly, but then very quickly. 
Even for very high noise levels w~ 1 — (100^iA") 2 rad -6 one 
may obtain a 3-sigma detection with a total sky coverage 
of about (4 deg) 2 (which corresponds to only 7 independent 
fields) For relatively low-noise w~ - 1 = (2(ktA") 2 rad- 8 , on the 
other hand, we may obtain an 8-sigma detection with a sky 
coverage of around (10 deg) 2 . Moreover, in the low-noise 
case, detections of several sigmas are possible for very small 
sky coverage. 



sqrti.sky coverage) (degrees) 



Figure 16. Contours of the detection function £ for stCDM for 
an interferometer with a primary beam with 9 W = 2° , for varying 
noise levels w~ 1 and sky coverage. 



opposed to E w 77L 2 /(4tt) for sCDM). Even in the pres- 
ence of noise, wherever the signal dominates, the detection is 
better for stCDM. However, in noise-dominated regions the 
behaviour of the detection function for stCDM and CDM is 
very different and the changes introduced follow a different 
logic for single-dish and interferometers. 

For single-dish experiments the signal-dominated re- 
gion is greatly reduced in stCDM. Much smaller beam- 
sizes 8b are now required for any meaningful detection. As 
shown in Figs. |l4|, and [li] one would now need 9b < 0.3° 
and 9 b < 0.25°, for noises w" 1 = (25^) 2 (deg) 2 and 
w^ 1 — (60/iif) 2 (deg) 2 respectively, in order to obtain any 
reasonable detection. Again one can plot an ideal scanning 
line in the beam/coverage sections defined by a fixed noise 
The ideal sky coverage is much smaller for stCDM 
than for sCDM. In general the countours of E for stCDM 
compared to sCDM are squashed to lower 9b, lower L, and 
achieve higher significance levels, with steeper slopes. Fol- 
lowing an ideal scanning line for any fixed w -1 one reaches 
a maximal detection allowed by the given level of noise which 
is always better for stCDM than for sCDM. This maximal 
detection is normally obtained with a small sky coverage, 
and infinite resolution. Nevertheless, one falls short of this 
maximum by only a few sigma if the resolution is about 
9 b = 4' - 6'. From Fig. |l|, for w" 1 = (25^iC) 2 (deg) 2 , one 
may now obtain a maximal 43-sigma detection for an ideal 
scanning area of L = 65 degrees. If 9b =0.1° a 36-sigma 
detection is still obtained. We also see that a beamsize of 
9b < 0.25° is required to obtain a 3-sigma detection (with 
L — 4 degrees), and a 10-sigma detection can be achieved 
only with 9t « 0.15° (and Li = 18 degrees). All-sky coverage 
for an experiment targeting stCDM is generally inadvisable, 
and it would only be optimal for the unrealistically low levels 
of noise vj- 1 < (llfiK) 2 (deg) 2 . 

For interferometers the signal dominated regions are 
similar for sCDM and stCDM (compare Fig. |l| with 
Fig. flq). The borderline detection regions (where noise dom- 
inates), however, have expanded. This is because for stCDM 
we may make the primary beam smaller without running 
into convolution problems, which in turn reduces the in- 
strumental noise. The high noise region is now w _1 > 
(150^) 2 rad -6 . There one should observe only one or two 



7 CONCLUSIONS 

We have studied the significance of secondary peak detec- 
tion for sCDM and stCDM in a large parameter space of 
experiments, including interferometers and single-dish tele- 
scopes, and have adopted a broad minded attitude towards 
sky coverage. If point source subtraction is to be done in par- 
allel with a CMBR experiment (so as to account for variable 
point-sources), howev er, a large sky-coverage may never be 
possible (see Section 4.3). This detail, often overlooked in 
satellite proposals, could then radically undermine a large 
number of estimates assisting experimental design. We will 
however not dwell on this awkward possibility, but while 
keeping it in mind, shall consider unrestricted sky coverage. 

The results obtained are reported in Sections ^ and ^ 
and stress the contradictions of an all-purpose experiment. 
If the low-/ plateau of the spectrum is a theoretical target 
then one needs all-sky coverage, and satellite single-dish ex- 
periments are to be favoured. As shown in Section |E| even if 
one wishes to study Doppler peak features for sCDM all-sky 
coverage might still be preferable. Depending on the noise 
levels, a large sky coverage might be desirable, even for a 
resolution of about 9b = 0.4 — 0.5°. 

Our work shows how such a design relies heavily on the 
assumption that the signal is in the vicinity of sCDM. If 
instead one is to test the high-i opposition between low 
CDM and cosmic strings, then we have seen that single-dish 
experiments are required to have rather high resolutions. 
Interferometers appear to be less constrained, providing 2- 
3 sigma detections under very unassuming conditions, with 
rapid improvements following further improvement in ex- 
perimental conditions. Furthermore, in this context, all-sky 
scanning is not only unnecessary, but in fact undesirable. 
The best scanning is normally achieved with deep small 
patches. These two features contradict sharply the ideal 
experimental design motivated by the standard theoretical 
gospel. 

Overall the parameter space of successful stCDM de- 
tections seems to increase for interferometers and shrink for 
single-dish when compared with sCDM. We believe that a 
variety of contrasting experimental techniques may equally 
well find their niche as regards important theoretical impli- 
cations. 
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